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Abstract. 

We obtain in exact arithmetic the order 24 linear differential operator L24 
and right hand side of the inhomogeneous equation £24 (3?^) = B' 6 ', 

where <E>( 5 ' = x' 5 ' — X^ 3 '/2 + 3C /120 is a linear combination of n-particle 
contributions to the susceptibility of the square lattice Ising model. In Bostan 
et al (J. Phys. A: Math. Theor. 42, 275209 (2009)) the operator L 24 (modulo 

a prime) was shown to factorize into ■ L^ 5 ^ \ here we prove that no 

further factorization of the order 12 operator k'o"' is possible. We use the exact 
ODE to obtain the behaviour of x' 5 ' at the ferromagnetic critical point and to 
obtain a limited number of analytic continuations of x^ 5 ' beyond the principal 
disk defined by its high temperature series. Contrary to a speculation in Boukraa 
et al (J. Phys. A: Math. Theor. 41 455202 (2008)), we find that x (5) is singular 
at w = 1/2 on an infinite number of branches. 
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1. Introduction 

The story of the zero-field susceptibility x °f the two-dimensional Ising model is a 
landmark saga of mathematical physics. A recent review of the highlights can be 
found in pQ. While a closed form expression for the susceptibility still eludes us, we 
possess an enormous amount of exact or extremely precise numerical information. This 
largely derives from two complementary approaches. One approach involves studying 
the series expansion of the susceptibility. Since the work of Orrick et al [2] we have had 
available a polynomial time algorithm, now of complexity 0(iV 4 ) for a series of length 
N terms. From the point of view of an algebraic combinatorialist, this comprises a 
solution, and many questions about the asymptotics, and about the scaling functions 
have been answered by analysis of the very long series we now have available-currently 
more than two thousand terms in length. 

The difficulty in proceeding further with this approach is that we have no idea 
what the underlying closed-form solution looks like, except that it is known, or more 
precisely universally believed, to be non-holonomic. The alternative approach is to 
express the susceptibility as a form-factor expansion. This approach was initiated 
more than 30 years ago by Wu et al [3J. In this representation, the susceptibility is 
written 

^x = (i-t) 1/4 Ex (2n+1) (1) 

n>l 

for T > T c , where t = sinh 4 ( J/k B T). For T < T c a similar expression with even 
superscripts prevails. The advantage of the form-factor approach is that each term in 
the sum is holonomic. This means that, with sufficient computational resources, and 
sufficient ingenuity, each term can be found. To date, the first six terms have been 
found, in the sense that their defining ODE has been obtained, either totally, or modulo 
a prime. We also have precise integral representations for the form-factor terms, and 
a Landau analysis of the integrands can provide information as to the distribution of 
singularities in the complex plane. Indeed, it was just such a study by Nickel [H [5] 
that gave convincing evidence of a natural boundary in the total susceptibility, thus 
supporting an earlier but weaker argument of Guttmann and Enting 6 that the total 
susceptibility was non-holonomic. 

While an exact solution for the Ising model susceptibility may be impossible and 
is certainly beyond reach at present, one might hope to obtain a complete picture 
of the singularities of the susceptibility. Indeed, this more limited goal has been the 
main motivation of the recent studies of the individual form-factor terms. For this 
the ODE and Landau analysis approaches are complementary; the Landau analysis 
provides necessary but not sufficient conditions J\ while the ODE, even if only in 
mod prime representation, can show which Landau singularities are to be excluded. 
The most detailed study in this regard is that of the five particle contribution x^ 
to the susceptibility initiated by Boukraa et al [8] and followed by Bostan et al [9]. 
The present paper is an attempt to address a number of issues left unresolved in these 
papers. 

A brief summary of the parts of [5] and [5] relevant here is as follows. In Boukraa et 
al [5] series in w modulo a prime to 10000 terms were given and shown to be adequate 
to find the order 33 Fuchsian differential equation^], modulo a prime, L^x^) = 0. 

X The notation used here is that in [5] : variables w and s are useful for high temperature expansions 
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Subsequently in Bostan et al [5] , the complexity of this ODE was shown to be reducible 
to an inhomogeneous equation 

L 24 ($ (5) ) = E<* 5 \ (2) 
where "J^ 5 -* is the linear combination of 5, 3 and 1-particle contributions 

$(5) = £(6) _£(3)/ 2 + x (1) /120. (3) 
The right hand side of ([2]) which satisfies L^(E < - 5 ' > ) = is of the form 

E {5) = w ■ [(1 - 16 w 2 ) 3 ■ P 4 ,o ■ K 4 + (1 - 16 w 2 ) 2 ■ P 3>1 ■ K 3 E 
+ (1 - 16w 2 ) ■ P 2 , 2 -K 2 -E 2 + Pi, 3 • K ■ E 3 
+ P ,4 • E 4 ]/ (1 +Awf/ '(1- 16 w 2 ) K , (4) 

where if = K(Aw) and i? = £(4w) are complete elliptic integrals and the 
Pi j = Pij(w) are polynomials. The degree of the polynomials and the denominator 
power k in ((4]) depend on the representation of L24 in ([2]). In the case that L24 
is minimum order 24, L24 is of degree 888, while the Pij are then of degree (at 
most) 904 with k = 8. In [9] results were only reported for a non-minimum order 
representation modulo a prime. Furthermore it was shown that L24 could be factored 
into L^f^ with L^' 8 being reducible into several smaller factors, all but 

one known in exact arithmetic. The question of whether could be factored was 

left unresolved. 

The results reported in [8] and [9] are an impressive example of what can be 
obtained by calculation modulo a prime. However there are limitations. Knowing L24 
and E^> modulo a prime enables one to deduce only possible local singularities of 
X^ and not its global behaviour. For example, one cannot determine the amplitudes 
of the singularities in x^ at the ferromagnetic point at s = 1. The leading and first 
correction term amplitudes were estimated in [8] but this was based on an analysis 
of the 2000 term exact integer series that had been obtained from multiple modulo 
prime series by the Chinese remainder theorem. 

In the present paper we determine the minimum order L24 and E^ in exact 
arithmetic. The results can be found on the website [10] . We have also constructed 
the minimum order exact L29 defined by i^g^ 5 ') = from these results. We have 
not actually needed this operator but provide it nevertheless for those who might find 
it of interest and in addition we also report in [TU] the exact integer coefficients of 
to 8000 terms that we generated directljH from @. The global information provided 
by @ enables us to prove that L^f^ defined by i 2 4 = L^f^ i[ r 2 Isht ' cannot be 
factored. We confirm the 500 digit amplitude of the leading ferromagnetic singularity 
of reported by Bailey et al [11] and also some |s| = 1 circle singularity amplitudes 
derived by Nickel [H [5]. 

An important question left unresolved in [8] was whether the singularity of the 
ODE for x^ at the zero of the head polynomial at w = 1/2 was in fact a singularity 
of x^ on some branch of the function. Appendix D in [5] provided an example for 
which not all the zeros of the head polynomial of the ODEp satisfied by an integral 
were Landau singularities of the integral. However, the integral for x^ wa s deemed 

with w = s/2/(l + s 2 ). The high temperature ^( 2n + 1 ) and ^( 2n+1 ) are related by s ■ ^-( 2n + 1 ) = 
(1 — s 4 ) 1 / 4 ■ ^( 2n + 1 ) so that ^( 2n + 1 ) has the "simpler" divergence oc 1/(1 - 4«i) oc 1/(1 — s) 2 at the 
ferromagnetic critical point. 

§ For the required x' 3 ) see |Appcndix B] and references therein. 

t We explicitly exclude zeros associated with apparent singularities. 
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too complicated in [5] to perform a complete Landau singularity analysis leaving only 
the conclusion that w = 1/2 was very likely not a Landau singularity of x^ ■ Here 
we perform some analytic continuation of x beyond the principal disk \s\ < 1 and 
onto other branches. Our exploration, while not exhaustive, is complete enough to 
show x^ has the singular behaviour (1 — 2w) 7 ^ 2 with non-vanishing amplitude on 
an infinite number of branches. We have not made any progress in identifying the 
Landau integrand singularities that give rise to this behaviour. 

In section [2] we describe briefly how we can combine results reported in [9] with 
multiple modulo prime <fr( 5 ) series of length at most 4000 terms to obtain the exact 
L24 and E^ in ([2]). This allows us to determine the series expansion of x*- 5 ' a t the 
ferromagnetic point by numerical matching of solutions. The results are reported in 
|Appcndix B| Section [3] is the outline of the proof that defined by L24 — 

L^2 ft ' • L^2 &ht ^ cannot be factored. Finally in section @] we describe the analytic 
continuation of <J>( 5 ' we have performed to obtain information on the behaviour of 
£(5) a t w = 1/2. 

2. The ODE for $( 5) in exact arithmetic 

We are looking for the unique minimum order L24 and associated E^ satisfying @ 
for a given $' 5 '. This is to be done modulo a prime for enough prime^j] that the exact 
integer L24 and E^ can be reconstructed by the Chinese remainder theorem. To 
generate the required $^ 5 - ) in (j3|) directly from the defining x^ integrals is impractical 
as the minimum order L24 is of degree 888 implying 25 x 889 = 22225 unknown 
coefficients in £24- Added to this are the five polynomials of degree 904 in E^- 5 ', i.e. 
another 5 x 905 = 4525 terms. Finding all coefficients is a straightforward linear 
algebra problem but requires that we have each $( 5 ) modulo prime series to about 
26800 termfl 

The operator L defined by L((f>) — for a given <j) is not unique if one does not 
require L to be of minimum order. The advantage of seeking a non-minimal order L 
is that the number of unknown coefficients to be found can drop dramatically. For 
example, only about 6200 terms are needed to obtain the non-minimal L29 of order 
51 and degree 118 used in [9]. The analogous effect occurs for the inhomogeneous 
equation ([2]) and we find that $^ 5 - ) series of length only about 5300 terms are needed 
when L24 is chosen of order 42 and degree 103 leading to polynomials of degree 155 
m E^ in (0} with now k — 26. Note however, as observed in [8], that the integer 
coefficients in non-minimal order operators can be outrageously large and we expect 
that Chinese remainder reconstruction of the non-minimal order L24 would be nearly 
hopeless. Instead, the utility of a modulo prime non-minimal L 2 4 hes in its use as a 
recursion device to extend directly generated (short) $^ 5 ' ) series to series of sufficient 
length, i.e. 26800 terms, so that the minimum order i 2 4 ar >d associated E^ can be 
found. Such extension requires completely negligible computer resources. 

A useful variant of the above approach is to use the non-minimal L24 found as 
described above to generate a long series for a solution ^24 satisfying L24(5 , 24) = 0. 
This series need only be about 22300 terms, long enough to enable reconstruction of 

% We have found that about 90 primes p < 2 15 are sufficient. 

X An alternative to Q is L 2 9(&^) = but this requires $( 5 ) series to 30 X f238 = 37140 terms. 
The change of our ODE problem from homogeneous to inhomogeneous form is a substantial reduction 
in complexity. 
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the minimum order L 2 4- Finding the coefficients in E^ from £ 2 4( < I >< ' 5 )) is then a 
separate and simpler problem. 

A further reduction in the length of the <fr( 5 - ) series to be directly generated can be 
obtained if one knows a factorization of L24 with the right division operator in exact 
arithmetic. This situation can be realized given the four modulo prime series reported 
in [5]. We have i 2 4 = £12 ' £ 12 sht ' 1 ■ Knowing in exact arithmetic allows 

one to obtain any non-minimal order representation of L^' 8 modulo any prim jjjl 
and then W = L^ ight) ($( 5 )) modulo a prime from the directly generated $( 5 ). If we 
choose our representation of the known as order 18 and degree 42, and the 

unknown as order 32 and degree 89, then the polynomials in = Lf 2 eft) (*) 

are of degree 199. This gives 33 x 90 + 5 x 200 = 3970 unknown coefficients in E^ 
and to be determined implying that 4000 terms of the directly generated x^ 

series is more than adequate. Since the series generation of x described in [8] is an 
0(N 4 ■ ln(iV)) process, this represents a 3-fold reduction in computer time from that 
needed using an unfactored L24 or a 6-fold reduction relative to the unfactored -L29 
approach. 

We have generated the x^ 5 "* series to a minimum 0(w 4000 ) for 90 primes p < 2 15 
and from these generated L 2 4 modulo a prime as described above. To obtain the 
exact i 2 4 is then a problem of rational reconstruction but it is relatively easy from 
just a few terms to guess a normalization factofl that converts the problem to 
integer reconstruction by the Chinese remainder theorem. The integer coefficients 
in L 2 4 are observed to typically have very large powers of 2 as factors which one 
can determine by a process of trial division by 2 k . If k is chosen too small, 
Chinese remainder reconstruction with a fixed number of primes might fail because the 
unknown coefficient is too large while if k is chosen too large there is failure because 
the coefficient is no longer an integer. With 84 primes we find an intermediate k range 
that yields a consistent integer reconstruction for every coefficient in L 2 4- With 90 
primes we have a large number of consistency checks that leave no doubt that our 
reconstruction is exact. We have also confirmed the apparent singularity constraint 
equations (A. 8) in [8] are satisfied by our reconstructed L24 in all cases, that is, 19849 
satisfied conditions on 22202 non- vanishing coefficients in i 2 4. 

In view of the "massive" calculations required to find L 2 4 in exact arithmetic 
it is natural to ask for some further mathematical and numerical checks of the 
correctness of the operators L 2 4 and £12 • First of all we have checked directly 
that does indeed right divide L 2 4 in exact arithmetic. Secondly, we have 

confirmed that the exponents of the operators L 2 4 and are rational numbers 

and in agreement with our previous massive numerical calculations 8 . Thirdly, the 
minimal order operator L 2 g corresponds to an integral of an algebraic integrand, and 
it is therefore, as mathematicians say a Period ("Derived From Geometry" [12]): L 2 g 
is thus, necessarily, globally nilpotent. This is a stronger constraint than being a 
Fuchsian operator (with integer coefficients) having only rational exponents. Since 

ff This is most easily done by using the minimum order L^ 5 ^ to generate a series solution S12 
1/^i g (S12} = but does not satisfy L n (S\2) = for any n < 12. The non-minimal 
- / I2 lgh *^ ' s then found from the series using, for example, the matrix code described in section 3 
of 0. 

f This becomes the value of the head polynomial at w = and is 2 12 -3 13 -5 7 -7 6 ■ ll 4 -23-29-7225564279 
= 4235287273136998077435560752320000000. 
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L 2 4 is a factor of L 2 g it must also be globally nilpotent, and likewise the left factor 
L^^ must be globally nilpotent. We have verified that L 2 4, and £i 2 Isht '[l 

are consistent with globally nilpotent operators. This is a very strong indication that 
our exact expressions for L24 and are correct. To check global nilpotence 

numerically requires one to calculate the p-curvature and check that it is zero for 
almost all primes. In practice one can obviously only do this for the first few primes. 
The primes used for Li 2 were those smaller than 30 while primes less than 10 were 
used for L 2 4- 

2.1. Computational details 

As shown in [5] the calculation of a series for x^ is a problem with computational 
complexity 0(N 4 In N). In [5] we initially calculated x^ to 10000 terms which 
required some 17000 CPU hours on an SGI Altrix cluster with 1.6 GHz Itanium2 
processors. From this single series we could already exactly identify a simple right 
divisor of L29 , and using this factor we were able to find a solution modulo a second 
prime using a series of 'just' 5600 terms. Expanding the series to order 5600 took 
around 1560 CPU hours on the Altrix cluster. The series modulo these two primes 
then sufficed to find a larger right divisor of L 2 g in exact arithmetic, and using this 
factor we found that only 4800 terms would be required to find solutions for any 
subsequent primes. 

Shortly after these developments a new system was installed by the National 
Computational Infrastructrure (NCI) whose National Facility provides the national 
peak computing facility for Australian researchers. This new system is an SGI XE 
cluster using quad-core 3.0GHz Intel Harpertown cpus. Our code runs almost twice 
as fast on this facility compared to the Altrix cluster. We then used this system 
to calculate the series for %( 5 ' to order 4800 (this took about 450 CPU hours) for 
a third prime, which again allowed us to find an even larger right divisor of L 2g in 
exact arithmetic. This larger operator reduced the required number of terms to 4600 
and we then calculated a series for a fourth prime to this order using some 380 CPU 
hours. These calculations gave us results for 4 different primes and allowed us to 
reconstruct the factor L^ 2 leht ^ in exact arithmetic as begun in [5] and completed here 
in |Appcndix A| 

It was only after this that we realised that the inhomogenous equation ^ could 
be used, as detailed above, to obtain simultaneous solutions for L 2 4 and E^ using as 
few as 4000 terms. We then calculated series for x^ to order 4000 for a further 86 
primes with each prime requiring about 215 CPU hours. 

The above timings make it clear that the reconstruction of L24 and E^ is 
a computationally expensive project. The main computational effort is the direct 
calculation of the series x^ modulo the required number of primes. It is therefore 
of some practical interest to estimate the number of primes required for the exact 
reconstruction based only on some partial reconstruction. We focus here on the 
coefficients of the head-polynomial of L 24 and denote by c„ the n'th coefficient 
(n = 0, . . . , 888) after stripping it of any factors of 2 as mentioned above. Then 
r n = ln(c„)/ln(30000) is a rough measure of the number of primes needed to 
reconstruct c n . From our previous reconstruction of L^ 5 ^ we noticed that the 
corresponding r n are given roughly by a quadratic function of n. Thus one can estimate 

X However, this operator is automatically globally nilpotent courtesy of its direct sum construction. 



ODE in exact arithmetic 7 



100 




300 600 900 

n 

Figure 1. Estimates r n = ln(c„)/ ln(30000) for the number of primes required 
to reconstruct the coefficients c n . 

the number of primes needed to reconstruct the full L24 from a partial reconstruction 
since r n for n near or 888 can be obtained using many fewer primes. In Fig. [T]thc 
lower 'curve' is the actual data from the coefficients of the head-polynomial of L24, 
while the upper solid curves are quadratic fits to the data based on the part data 
set, from top to bottom, r„ < 30, 40, 50, and 60. The pertinent point being that 
after doing the calculation for some 30 primes we estimated that the reconstruction 
was likely to succeed with no more than 100 primes and was therefore achievable in 
practice. 

With L24 and E^> known in exact arithmetic it is easy to use ^ to calculate the 
exact series for <J>^ 5 ^ to high order. Specifically our reconstruction means that we know 
the coefficients a,ij of the polynomials in the operator L24 and the coefficients of the 
polynomials P^j of (j4]) exactly. The latter allows us to easily calculate the coefficients 
e„ of E^ by using the (simple) recursive formulae for the elliptic integrals E and K. 
From © we have explicitly, by equating the coefficients of i", that (recall that L24 is 
expressed in terms of the differential operator x-^) 

M D 

^ ^ a-ij (n - j) 1 c n -j = Acn + B = e„, (5) 

i=0 j=0 

where A and B are integers (depending on n). The coefficients c„ of can thus 
be calculated recursively and from ([3]) we can calculate the coefficients of (with 
the coefficients of calculated using the ODE from [13]). We have calculated the 
coefficients of x' 5 ^ up to order 8000 and they can be found in jlOj . 

Finally we decided to calculate the minimal order operator L29 explicitly; it is 
given implicitly by ([2]). This can be done in a variety of ways. The most obvious way 
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is to use ([5]) to extend the series for <E>( 5 ' to high enough order (30 x 1238 = 37140, 
since the minimal order L 2 g has degree 1237) and then use the matrix code of [5] 
to calculate the ODE corresponding to L 2 g modulo a sufficient number of primes to 
reconstruct L29. However, computationally it is easier to first calculate modulo a prime 
the minimal operator L5 annihilating E^ and then form the product L5 • L24 modulo 
a prime. The minimal order operator for L5 has degree 4489 so 6 x 4490 = 26490 
terms of E^ is required. Obviously, since the minimal order L24 has degree 888, the 
product L5 -^24 has degree 5377 meaning that there is a common factor of degree 4140, 
which we must discard in order to calculate the degree 1237 polynomials of L 2 g. The 
minimal order L§ was calculated (for each prime) using the matrix code of [8]. The 
product L5 • L24 was then calculated modulo a prime using Maple and the common 
factor can then be divided out modulo a prime. The bottle-neck in this calculation is 
the use of the matrix code of [5] which has computational complexity 0(N 3 ). It was 
for this reason that we chose the "indirect" route of going through L5 to get I/ 2 g. 

3. Proof that d° es n °t factorize 

The factorization of L 2 g defined by L29(^ 5 ^) = described in [9] relied mostly 
on the testing of series S(w) generated by L2g(S) = modulo a prime around 
w = 0. If a particular series is annihilated by an L n for order n < 29 then L n 
right divides £29- Depending on the singularity exponents, a series solution^] S(w) = 
w q ■ (1 +oti w +ct2 w 2 + ■ ■ •) with fixed q, might be uniquely determined by L2g(S) = 
or might contain one or more arbitrary rational coefficients on. In the latter case the 
series may be a generator for a right division operator only for a particular choice 
of constants on and the problem then is how these particular values might be found. 
If there is only one arbitrary coefficient a in the series S(w) an exhaustive search is 
possible in a modulo prime, p, calculation because one need then only test p separate 
series with a an integer satisfying < a < p. If there is more than one arbitrary cti in 
the series S(w) such brute force "guessing" is no longer practical due to computational 
time constraints. It is this that prevented the authors of [9] from deciding whether 
is factorizable. Let us also note that we were not able to perform a straight 
formal calculation factorization of L^ff using its expression in exact arithmetic and 
our attempted calculations failed on a computer with 48 Gb of memorjQ. 

As pointed out in [9] there is nothing special about the singular point w = 0. 
And indeed, a series solution about a singular point w s ^ might also lead to a right 
division operator. For example, the unique singular series 

x~ 7/4 ■ (1 + 387 x/80 + 72103 x 2 /23040 + 2054561 x 3 /1597440 + 0{x 4 )) (6) 

with x = w — 1/4 is annihilated by an order three operator. With x = w — 1/2, the 
unique singular series 

x 7/2 ■ (1 -Alx/6 + 26557 x 2 /792 - 8692015 x 3 /61776 + 0(x 4 )) (7) 

yields an order 6 right division operator Lq. By proceeding through a sequence of 
such solutions, one can eliminate much if not all of the p-fold searching described in [9] 

f If the solution contains powers of ln(w) the S(w) here is to be interpreted as the coefficient of the 
highest power of ln(«>). 

If More precisely the Maple 13 command DFactor(*,onestep) was not able to yield a conclusive 
answer for such complicated large order operators. Seeking for a left division of we had 

similar inconclusive Maple calculations on ad^Lj 1 ^'), the adjoint operator of L^ lt \ 
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to achieve the factorization L 29 = L 5 ■ L 24 = L 5 ■ ■ L^" sht ^ and the additional 

factorization of Lfj r 2 sht ' ) . 

The points chosen for series expansion need not be restricted to the rational head 
polynomial roots as in the equations ^ and ([7]) above. The factor 1 + 3w +4w 2 , 
which appeared already in the head polynomial of the L-j that annihilated x*- 3 ', has 
"accidental" modulo prime factorizations for roughly half the primes close to 2 15 . For 
example, with prime p = 32719, one has 1 + 3w + Aw 2 — A (w — 8973) (w — 31925) 
modulo p. We can write x = w — w p with w p either 8973 or 31925 and obtain 
solutions about x = satisfying L2i(S(x) -In 2 (a;) + R(x) •ln(x) + Q(x)) = modulo 
p, where R and Q are regular at x — and the series S(x) — x + 0{x 2 ) is unique. 
Testing shows that S(x) is annihilated (modulo p) by an order 3 operator L^{x). As 
we show in general in | Appendix C| one can obtain from L^(x) a right division (modulo 
p) operator L^(w) and then from multiple modulo prime calculations, an exact right 
division L^(w) by the Chinese remainder theorem. This reconstructed £3(wj3 has the 
factor 1 + 3w + 4w 2 in its head polynomial in spite of the fact that w p is clearly not 
one of the roots (—3 ± i-\/7)/8 modulo p, of 1 + 3w + Aw 2 . 

Testing at points other than w — also enables one to exclude certain series 
solutions as generators of right division operators. For example, the case 3 polynomial 
1 — 7 w + 5 w 2 — A w 3 which is a factor of the head polynomial of 1/24 has the modulo 
prime, p = 32749, factorization 32745 • {w 2 + 11821 w + 10836)(w - 3635). Define 
x = w — 3635. Then the singular solution modulo po about x = is S(x) hi(x) +R(x) 
with 

S(x) = x 5 + 13877x 6 + 9339 a; 7 + 25021a; 8 + 21884 a; 9 + 0{x 10 ) (8) 

unique and R(x) regular at x — 0. Testing the series ([8]) shows there is no L n (S) = 
modulo po for any n < 24. Thus there is no operator of order less than 24 that has 
x = w — 3635 as a factor (modulo po) of the head polynomial and more generally, 
1 — 7 w + 5 w 2 — 4 w 3 as a factor. This also implies that any solution S that is singular 
at some root of 1 — 7w + 5w 2 — Aw 3 and satisfying L^^{S) = 0, cannot also be 
a solution of an operator of order less than 12 that right divides L^^. The same 
conclusion is reached for the remaining case 3 and both case 4 polynomials for x' 5 ' 
from Appendix C in [5]. 

If we only knew i 2 4 or modulo prime for a few primes, the above 

information about singular solutions associated with case 3 and 4 polynomials would 
not be particularly useful for finding or excluding factorization. However, with the 
exact i 2 4 available one has global information and can match series solutions about 
w = to solutions about other singular points of L 2 4. In particular, if one can 
show that every series solution S(w) about w = satisfying L^^iS) — is 
singular at some root of the case 3 or 4 polynomials then one has proved that 
does not factorize. Our demonstration that this is the case uses the singular point 
w s = 0.15853 • • • which is a root of 1 — 7w + 5w 2 —Aw 3 . This is a particularly 
convenient point as it is the closest root of the head polynomial of L 2 4 to both w = 
and w — 1/4. 

\ It is equivalent to the product Z2 ■ N\ which was shown on general grounds in [9] to right divide 
£29- 

f We follow the notation of Appendix C in [8] . 

|| Equivalently every solution of L24 that is not a solution of L^" 8 * 1 *'- 
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We begin the demonstration by studying the two linearly independent solutions 
of the form Si = Ai(w) ■ ln 3 (u>) plus terms with lower powers of hx(w). The two series 

A 1 = 9 w + 261 w 3 + 1845 w 4 + 7046 w 5 + 42771 w 6 + 145980 w 7 

+ 785528 w 8 + 2536628 w 9 + 12800309 w 10 + 38627228 w 11 

+ 187738058 w 12 + ■ ■ ■ + a hn ■ w n H , (9) 

A 2 = 27 w 2 +102w 3 + 270w 4 +2164w 5 + 5532w 6 + 43722 w 7 
+ 132130w 8 + 922108 w 9 + 3158590 w 10 + 19690882 u; 11 
+ 72977164 w 12 + • • • + a 2 , n ■ w n + ■ ■ ■ , (10) 

satisfy Li^yAi) = but not (Ai) = and a ratio test shows they have radii 

of convergence \w\ = w s = 0.15853 • • •. Thus both A\ and A2 are singular at w = w s 
and cannot be generators of an L n , n < 24, that right divides £24- Whether a linear 
combination of A\ and A2 leads to a right division operator is now determined as 
follows. 

Near x = where x = w s — w the Ai must be of the form Ai = Bi ■ 
f(x) ■ hi(x) + gt(x) where f(x) and the gi(x) are all regular at x = 0. The series 
A(w) oc £?2 Ai(w) — B\ A2(w) will be regular at w = w s and if the amplitude ratio 
B1/B2 is rational, then A(w) might be a candidate generator for a right division 
operator. In principle, the amplitudes Bi could be found by matching series solutions 
about x — to those about w = but since we only need the B1/B2 ratio, a 
simpler procedure that utilizes only the coefficients in (|9I10|) is possible. We note that 
of the remaining singularities of Ai, the nearest to w = are at \w\ = 1/4. This 
implies that B1/B2 — cti, n /oi2,m a ratio of coefficients from (|9ll0p . to an accuracy of 
order (Aw s ) n « 0.634". Thus the problem reduces to searching for a (small) rational 
B1/B2 from a sequence that converges exponentially. This can be done by expressing 
ct\,njoL2,n as a continued fraction. We observe that, for n greater than some fixed 
no, a particular term in the continued fraction grows exponentially which is a clear 
indication that in the limit n — > 00 the continued fraction terminates and is the 
(small) rational B1/B2 = —637/228. Our result for the linear combination series is 
then 

A{w) = 2052w + 17199w 2 + 124482u; 3 + 592650w 4 + 2984956w 5 + 0{w 6 ) (11) 

which is confirmed to have a radius of convergence |w| = 1/4. However, testing 
the series (|11[) shows it is not annihilated by any L n , n < 24, and thus we have 
eliminated the only possible linear combination candidate series for a right division 
operator. 

While direct testing is an easy way to exclude A(w) in (JTTJ), such direct testing 
is impractical as a general method for excluding the many possible series that arise 
in the remaining part of our proof. Instead we supplement direct testing by a method 
that relies on the global information provided by the exact £24- As an illustration 
of this method, consider again A(w). If one matches ([TT|) to series about w — 1/4 
one finds that A(w) contains, as the leading logarithmic function, A s ■ ln 2 (y) where 
2 y = 1 — 4u> and 

As = 1 - 21469 ?//640 - 1489293 ?/ 2 /81920 + 229328363 y 3 /i0485760 + o(y 4 ) (12) 

A ratio test on the coefficients in (|T2"j) shows A s has radius of convergence \y\ = 
1/2 — 2 w s and thus is singular at w = w s . This in turn implies that there are 

§ This was the problem encountered in [9]. 
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branches of the function A(w) on which it is singulaiQ at w — w s . By forming the 
linear combination A(w) oc B 2 Ai{w) — BiA 2 (w) we only succeeded in forcing an 
"accidental" cancellation of the singularity at w = w s on the principal branch of 
the function. The singularity remains on at least some other branches and thus A is 
excluded as a generator of a right division operator for exactly the same reason as Ai 
and A 2 . 

The above argument has excluded, as generators of right division operators, eight 
of twelve linearly independent solutions satisfying (S) = 0. These we take to be 

Ci oc 42 ight) (^) with A from CP and C 2 oc L { ^ eht) (A 2 ) with A 2 from (10]) plus 
the six series with leading ln(ui) dependence C\ ■ \n p (w), p = 1, 2 and 3. Explicitly, 

Ci{w) = w 6 - 444W 7 / 11 + 275773109 w 8 /129360 + 19252320091 w 9 /194040 

- 964738631897 w 10 /388080 - 2082457681309 w 11 /27720 

+ 17517580633073581 w 12 /17075520 + 0( w 13 ), (13) 

C 2 {w) = w 6 + 403206 w 7 /1661 - 13446782071 w 8 /19533360 

- 2413114741889 w 9 /29300040 - 4359267083039 iu 10 /1065456 

- 875856906689449 w n /4185720 

+ 23619065101886078533 w 12 /2578403520 + 0{w 13 ). (14) 

Because L^ 2 &Ylt ^ does not have the factor w — w s in its head polynomial, the C(w) 
functions carry the same w = w s singularities as the A(w) from which they have been 
generated. Thus even though the series Ci(w) has radius of convergence |w| = 1/4, 
the analytically continued function C\(w) is still singular at w — w s on some other 
branches. The series C 2 {w) has radius of convergence \w\ — w s and is already singular 
at w = w s on the principal branch. 

It was shown in [9] that the remaining four solutions satisfying (S) — are 

of the form d(w) -ln(w) + Di(w) and Ci(u>), i = 3, 4, with the d and Di regular at 
w = 0. These two Ci, together with the two in (|13l[T4"|) . are linearly independent and 

can all be generated from the coefficient of ln(ui) in L^ 2 sht \s 2 4) where S24 contains 
four arbitrary constants and is of the form S24 = F(w) ■ In 2 (w) + G(w) ■ ln(w) + H(w) 
with F, G and H all regular at w = 0. We demand that F satisfied L^ ight) (f ) = 0. 
This guarantees that C(w) ■ ln(w) is the leading logarithm in Lj r 2 Isht ' ) ( l S'24) and 
simplifies the subsequent analysis. Only F and G are relevant for determining C 
and a possible choice is 

F = fti ■ (3177 w - 174840 w 4 - 817828 w 5 - 5829558 w 6 

- 25983762 w 7 - 142882882 w s - 620769318 w 9 - 3086072424 u; 10 

- 13199839762 w 11 - 62214586728 w 12 + 0{w 13 )) 
+ f3 2 - (3177 w 3 + 13803 w 4 + 74932 w 5 + 287997 w 6 

+ 1265280 w 7 + 4296418 w 8 + 17162736 w 9 + 48945231 w 10 

+ 173557768 w 11 + 284486847 w 12 + 0{w 13 )), (15) 

X These singularities can be reached, for example, by following a path along the real axis from w = 
to w just less than 1/4, circling w = 1/4 any number N of times, and then moving back to w = w a 
along the real axis. Since A s multiplies the leading (unique) logarithmic singularity at w = 1/4, 
there cannot be cancellation of the w = w B singularity of A s on all branches distinguished by N. 
§ For the explicit F in (115b . it happens that d_F/d/3i is annihilated by an L\\ and dF/d/32 by an 
Lg. This reduction from L12 plays no role in our subsequent arguments. 
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and 

G = Pi - (1604883673 w 8 /210 + 2823208099 w 9 /105 

+ 47115755881 w 10 /140 + 782148892459 w u /630 + 0{ w 13 )) 
- p 2 ■ (366106439 w 8 /1050 + 1576821038 w 9 /525 

+ 9206778909 w 10 /350 + 1049578781449 w n /6300 + 0( w 13 )) 
+ p 3 ■ (35 w 6 + 1223 w 8 + 1852 w 9 + 36064 w 10 

+ 96388 w 11 + 0(w 13 )) 
+ Pa • (105 w 7 + 304 w & + 3536 w 9 + 10192 u; 10 

+ 79089 w 11 + 0(w 13 )), (16) 

where the Pi are arbitrary constants. We will now show that no choice of these 
constants can yield a C = C(@i, p 2 , /?3, Pa), defined by 

4 r 2 ight) (F-ln 2 H + G ■ ln(iu) + H) = G-lm» + D, (17) 

that is a generator for a right division operator of . The argument is essentially 

that given above for the exclusion of Gi and C 2 in (fT31 [T4"| . In fact the demonstration 
has already been partially completed since, in terms of C{P\, P 2l P3, Pa), C\ oc 
(7(0,0,3,4) and C* 2 cx C(0, 0, 453, 1744). Since F in (TT5]) satisfies L^ 2 ight) (F) = 
it is not singular at w — w s — 0.15853 neither is ght) (F • ln 2 (w)). Thus it 
suffices to investigate G and if every G is singular at w = w s then so is C defined by 
(|17l) and we have proved does not factorize. 

A ratio test on the series coefficients in (|16l) shows the generic G has radius of 
convergence \w\ — w s and thus is singular at w = w s . But by the same analysis 
that led from the series (J9J [TO)) to the linear combination (flTT) , we can construct three 
G = G(pi, P2, P3 7 Pa) each of whose radius of convergence is \w\ — 1/4. These are 

G(105,0, 0,-1182781), G(0, 525, 0, -1443727), G(0, 0,3,4). (18) 

To the remaining linearly independent G — G(0, 0, 0, 1) one can add any combination 
of the three in ([T8"l) but this will not change its radius of convergence from \w\ = w s 
and remove the singularity at w = w s . In this sense G(0, 0,0, 1) is equivalent to 
G(0, 0,453, 1744) which corresponds to C 2 via (|17l) and is excluded as a generator of 
any right division operator. 

To determine the behaviour of the three G functions in (TT5| in the vicinity of 
y = where 2y = 1 — 4w we match the serie tQ 2F-lm» + G in w about w = to 
solutions S satisfying L 2 4(S) — about y = 0. Since one can show 

that F (and F ■ hi(w)) can contain only the first power of ln(y) near y = 0. Any 
ln 2 (y) or ln 3 (y) we find in the matching S can only come from G in the combination 
solution 2 F ■ ln(w) + G. Our matching shows the leading logarithmic dependencies 
of the three G's in (|18l) are respectively 

(-3420025/8192/tt 2 ) A s ln 3 (y), (-2473625/16384/tt 2 ) A s ln 3 (?/), 
(-125/16384/tt 2 ) A s ln{y) (19) 

where A s is given by (fl2|) . A linear combination of the first two G's in ([T8"l) 
can be constructed to eliminate the leading A s ln 3 (y) shown in (fT9|) and we find 

f It follows from L2i(F ■ ln 2 (ui) + G ■ ln(to) + H) = by analytic continuation around the w = 
singularity that also L2i{2F ■ ln(to) + G) = 0. 
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the resultant G(494725, -6840050, 0, 13236968029) has A s ■ ln 2 {y) as the leading 
logarithmic singularity. Clearly this remainingsingularity can now be eliminated by 
forming a linear combination with the last G ir|3] (|18[) . In summary, we have generated 
the three G(/3i, fc, Psi AO combinations 

G(494725, -6840050, 15276842775, 33606091729), 

G(0, 525, 0,-1443727), G(0,0,3,4) (20) 

which are irreducible in the sense that no further superposition can eliminate the 
A s ■ ln 3 (y) and A s ■ ln 2 (y) singularities of the last two G's while the first G is unique 
in that it contains neither A s ■ ln 3 (y) nor A s ■ ln 2 (?/) . The last G in (f2"01 we have 
already identified as being associated with C\ . We associate the middle G in ([20]) with 
G3 oc G(0, 525, 0, —1443727) which cannot be a generator of a right division operator 

of for exactly the same reason as C\. The first G in (j2T)]) we associate with 

G 4 oc G(494725, -6840050, 15276842775, 33606091729) and test C 4 directly. We find 
there is no operator satisfying L n {C$) — with n < 12 and this completes our proof 
that does not factorize. The explicit new C,-, supplementing those in (I13I14P 

are 

G 3 (» = w 6 + 281575923 w 7 /34167793 
+ 48755202697119 w 8 /8371109285 
+ 788146152364265 w 9 /5022665571 

- 48321460210711729 w 10 /33484437140 

- 1046678480403963299 w n /4783491020 

- 4705373665277858926411 w 12 /3683288085400 + 0( w 13 ), (21) 

C A {w) = w 6 - 124536 w 7 /649 + 2840488261 w 8 /508816 
+ 39013193251 w 9 /254408 - 24532098411899 w 10 /9540300 

- 6082145734733 w n /68145 

- 179169570633725593 w 12 /314829900 + 0(w 13 ). (22) 
3.1. More on the structure of the differential operator 

Once the differential operator has been proved to be irreducible, one may wonder 

whether this high order differential operator can nevertheless be built from factors of 
lower order. A high order differential operator can be irreducible and still result from 
"operations" involving differential operators of lower order since it may be a symmetric 
power of a lower order differential operator or a symmetric product of two (or more) 
lower order differential operators. 

The symmetric n'th power of a differential operator L q is the differential operator 
whose corresponding ODE annihilates a generic linear combination of the q solutions 
of L q to the power n. The symmetric n'th power of a differential operator L q of order 
q has order ^|' . 

The symmetric product of the differential operators L qi and L q2 of orders qi and 
q2, respectively, is the differential operator whose ODE annihilates the product of a 



1 The surprise, at least for us, is that all of the necessary combinations can be formed with rational 
amplitudes. 
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generic linear combination of the q\ solutions of L qi and a generic linear combination 
of the <? 2 solutions of L q2 . This symmetric product is of ordeiQ qi ■ q 2 ■ 

We use the notation [w p ] to indicate a series that starts as w p (const. +•••). In 
[5] it was shown that the formal solutions of at w = follow this scheme: There 

are two sets of four solutions (k — 6, 7) 

[w k ] \n(w) 3 + [w 5 ] ln(w) 2 + [w] \n(w) + [w], 
[w k ] \n(w) 2 + [w 5 ] ln(w) + [w], 

[w k ] \n(w) + [w], and [w k ] (23) 

and two sets of two solutions (k = 8, 9) 

[w k ] ln(w) + [w], and [w k ], (24) 

[w k ] \n(w) + H and [w k ] (25) 

We denote by BLn a set of solutions such as (|23l) containing n + 1 solutions with a 
logarithmic solution of maximal degree n. For the scheme above, we thus have two 
BL3 blocks and two BL1 blocks, and in each block there is also a non-logarithmic 
solution. 

We first consider the possibility that is a symmetric power of an operator 

of lower order. It is straightforward to see that the only possibility is that could 
be a symmetric eleventh power of a differential operator of order two with one BL1 
block. This possibility is ruled out, since there is no BLW block in the solutions of 

r (left) 
L \2 ■ 

Next for the possibility that is a symmetric product of differential operators 

of lower order (we consider only the cases where the product has the maximal order). 
There are three cases to consider: The symmetric product of differential operators of 
orders two and six (configuration denoted 2-6), three and four (configuration 3-4), 
or two, two and three (configuration 2-2-3). 

The symmetric product of two differential operators L\ and Li containing the 
blocks BLn\ and BLn2, respectively, should contain in its solutions the block BLn 
with n = ni + n-z- Since the differential operator L^^ contains two BL3 blocks the 
case 2 ■ 2 ■ 3 is ruled out. 

Let us detail the compatibility of the case 3 • 4 at w — 0. With two BL3 blocks 
in L^^ the only possibility is that the order three differential must have one BL2 
block: 

Si ln(ui) 2 + Sn ln(w) + Sw, 
Si ln(w) + 5 20 , 

Si (26) 
and the order four differential operator must have two BL1 blocks: 
Ti ln(iu) + Tio, 

Ti, (27) 
Vi \n(w) + Via, 

Vi (28) 



\ The order q = qi ■ q2 is for the generic case. In general the order of the symmetric product is 

9i + 92 - l < q < qi ■ qi ■ 
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It is a simple calculation to form the product of a combination from the set BL2 with 
a combination of the solutions from the two sets BL1. One obtains: 

Si • Ti InH 3 + (Si • Tio + S n ■ T x ) ln(u-) 2 + (S i0 • T x + S u ■ T w ) ln(to) + S w ■ T 10 , 

S x ■ T x Mwf + S u ■ T x ln(to) + Sio ■ T x , (29) 

Si • Ti ln(to) + S 2 o • Tx, (30) 

Si • Tx (31) 

and a similar set of four solutions with Vs instead of T's. These eight solutions 
correspond to the two BL3 occurring for L X 2 at to = 0. One obtains also 

Si • Tx \n(w) 2 + (Sx ■ T 10 + Tx • S 20 ) ln{w) + S 20 • T 10 , (32) 
Si • Tx ln(w) + Si • Tio (33) 

and two other solutions where V's replace T's. Subtracting (f3"2"j) from (f2"9"|) and (133)) 
from (|30|) . one obtains 

(Si • Tio + Tx ■ S 20 - Tx ■ Sxx) ln(u.) + (S 20 • Tio - S 10 ■ Tx) , (34) 
Si • Tio — Tx ■ S 2 o- (35) 

Note that in a set of solutions such as BUI above, the series Sn depends on S20 and 
Si and can be expressed as 

Sn = aSi + 2S 2 o (36) 

The coefficient 2 is generic for any order three ODE and a is a constant that depends 
on the ODE at hand. Inserting this Sn in (f34|). we can arrange to have the non- 
logarithmic series be the same as the series in front of the log. This is then one of 
the BL1 blocks occurring in L^^K The second BL1 block is obtained by considering 
([Mf with V's instead of T's. 

We have thus shown that under the hypothesis that is a symmetric product 

of two factors, the scheme of solutions at w = is compatible with the configuration 
3 • 4. Similar calculations show that the configuration 2 • 6 is also compatible and in 
this case the order two operator has a BL1 block and the order six operator has two 
BL2 blocks. 

Next we must check for each configuration (3 • 4 and 2 • 6) whether or not the 
symmetric product is compatible with the scheme of solutions for at those other 

singularities containing enough logarithmic solutions. We recall that the scheme at 
w = 00 is the same as the scheme at w — 0. The scheme of solutions of L^ 1 ^ at the 
singularity w = 1/4 is one BL3 and four BL1 blocks. One sees immediately that the 
configuration 3 • 4 is ruled out. For any set of solutions that we attach to the order 
three and order four differential operators, we end up with either more than one BL3 
or at least one BL2. The configuration 2 • 6 is acceptable, since in this case there can 
be one BL1 block in the order two operator and one BL2 plus three BLO blocks in 
the order six operator. It remains to be seen whether or not the configuration 2 • 6 is 
compatible with the scheme of solutions at the point w = —1/4, which is two BL2, 
one BL1 and four BLO blocks. Our checks show that the configuration 2 • 6 is ruled 
out. 

In conclusion we have shown that the differential operator LY^' isn't a symmetric 
n'th power of a lower order operator nor is it a symmetric product of two (or more) 
operators of orders q x and g 2 under the hypothesis that the order of the symmetric 
product reaches its maximum value qx ■ q2 = 12. It should be noted that our 
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considerations regarding symmetric powers/products (unlike the question about the 

factorization of L^^) do not require knowledge about Lif^ in exact arithmetic. The 
block structure of the solutions can be obtained from the operator modulo a prime. 

4. Analytic continuation of <£>( 5 ) and its behaviour at w = 1/2 

Our analytic continuation of $ < - 5 - 1 is limited to paths that follow the real w axis on 
the intervals [0, 1/4], then [1/4, (3 - %/5)/2], and finally [(3 - V5)/2), 1/2]. We allow 
any number of half-integer turns about the ferromagnetic point w = 1/4, that is, 
rotation by any angle 6 — n ir with n odd. This is followed by any rotation 8 — m n, 
m odd, around w = (3 — v / 5)/2. This point is one of the s-plane circle singularities 
discussed by Nickel [H[5]. The point w — 1/2 also maps onto the |s| = 1 circle but it 
is not a singularity of x when it is approached on any path that does not leave the 
principal disk |s| < 1. In terms of our n,m paths in the w-plane, the combinations 
n = m = ±1 are such principal disk constrained paths. Every other n, m combination 
is a path that reaches w = 1/2 on another branch. 

Our starting point for the analytic continuation is the $*- 5 ' series expansion about 
w = 1/4 in |Appendix B| For < w < 1/4, 2y = 1 — Aw is positive real and the 
half- integer turns about w = 1/4 that bring one to w > 1/4 simply requires the 
replacements 

V ~ > -V, ln(y/4) — > ln(|»|/4) -inn, (37) 

in (|B.3[) with n odd to be understood. The new (|B.3[) series generated with the 
replacements (|3"T)) is now to be matched to series in z where z = 3 — \/5 — 2 w = 5/2 — 
y5 — y. Although direct matching is possible, considerable improvement in the utility 
of the y series results by first making an Euler transformation by the replacement 
V 2//(l ~ v)- This has the effect of bringing the w = (3 — v / 5)/2 ~ 0.382 singular 
point closer to w = 1/4 while moving w s w 0.1585 further away. 

To generate $*- 5 ' ) series in z we must first analytically continue the elliptic 
integrals in ((4]). The replacements required for w > 1/4, with the same n as in 
1(57]). are 

K{Aw) — ► u[K{u) + \nK{u% 

E{Aw) -K{Aw) — > [E(u) -K(u) -\nE{u')]/u (38) 

where u = I /(Aw) and u' = y/l — u 2 . The new elliptic integrals ([38]) are easily 
developed as series in z from their defining differential equations set up as recursion 
relations. The remaining step of finding a particular integral of ([2]) together with all 
homogeneous series solutions is also straightforward. The ODE numerical recursion 
in z is not as unstable as the recursion in y noted in | Appendix B| Here one loses only 
about a factor 10 in relative accuracy for each two orders in z. 

The matched series in z is of the form A(z) ■ ln(z) + B(z) where A and B are 
regular at z = 0. Since and are n °t singular at this point, we can identify 
the singularity in $( 5 ) with that in x^- We find the leading singular term i43 

[(1 + 10n 2 + 5n 4 )/16] • [(25/693)(5 - V5)(2 + V5) n /(2 22 ir 2 )} ■ z 11 ln(z). (39) 
The n = 1 singularity, when mapped to s-plane variables, agrees with the sum 
contribution x^i + X() 5 -i from equation (14) in 4 . The new result in (|3"9")l is the 
branch dependent multiplicity (1 + 10 n 2 + 5n 4 )/16. 

| This is based on floating point results to about 300 digit accuracy. 
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To most clearly identify a possible w = 1/2 singularity in the A(z) and B(z) 
series we make another Euler transformation with the replacement z — > z/(l + 6 z/5). 
This moves the known singularity at w = 1/4 to the new z = 5(V5-1)/16 w 0.386 
and the potential singularity at w = 1/2 to the new z = — 5(16 - 5-y/5)/131 « 
—0.184. There is another potential singularity at a complex root of the case 4 
polynomial 1 — w — 3w 2 + Aw 3 . This maps to the new \z\ 0.387. Since the 
singularity of interest is about factor 2.1 closer than the next nearest, any singularity at 
w = 1/2 will be observable in an iV term Euler transformed z series with corrections 
of order (2.1)-^. This factor (2.1)"^ is also the bound we can put on any it) = 1/2 
singularity amplitude if a ratio test of coefficients does not indicate a singularity at 
z rj —0.184. We find the absence of such a singularity for the A(z) series which we 
have generated to length N = 1100. Because the A(z) series multiplies In (z), this also 
indicates that any possible singularity at w = 1/2 will have an amplitude independent 
of the index m specifying the logarithmic branch of the z — (w — (3 — y/E)/2) 
singularity. 

We find that for n — ±1, the B(z) series is also not singular at z « —0.184. For 
other n values a singularity is clearly indicated and by a coefficient ratio analysis of 
different n series completely analogous to what was done for the A\ and A2 series in 
(|9|10p we find that the singularities at w — 1/2 have amplitudes proportional to the 
branch dependent multiplicity factor (n 2 — l) 2 . A more detailed analysis involving 
explicit fitting of the B{z) series coefficients yields the singularity amplitude which 
we have confirmed by direct matching of the x^ series in z to ODE solution series 
about w = 1/2. The result for the singular part of x a t w = 1/2 is 

X2 S = [(n 2 -l) 2 (2V6)/(315^)] • (1-2H 7/2 

x [1 +41(1 -2w)/12 + 26557(1 - 2w) 2 /3168 + •••], (40) 

where the amplitude has been verified to our numerical accuracy of about 250 digits. 
We have not identified a Landau singularity associated with (|4Tfl) . Finding this Landau 
singularity in the x. integrand remains as the major unsolved challenge of this paper. 

5. Conclusion 

We have completed the quest begun in [8] for the exact integer arithmetic ODE 
satisfied by x^ 5 "*- While most explicit results are far too extensive to be published 
here, a selection can be found on the website [TU] . These include L24 and E^> defined 
in ([2])-((4]), the explicit factorization L24 = ■ L^ ebt \ the operator L29 and the 

exact integer high temperature series for x^ to 8000 terms generated from ([2])-(|4]) 
and numerical coefficients to 800 digits for the x series at the ferromagnetic point 
to supplement (|B.4|) . 

We have used the exact ODE to resolve at least one issue that was left undecided 
in [5]. In particular, we have shown in section [3] that the operator cannot be 

factored. The techniques we have described, both for factorization and proving the 
converse, are not all well known and we expect they will find application in other 
problems. In addition we have shown that isn't a symmetric n'th power of a 

lower order operator nor is it a symmetric product of two (or more) operators if the 
order of the symmetric product reaches its maximum value. 

An issue raised in [8] has only been partially resolved. No obvious candidate 
for a Landau singularity in the integral representation of x' 5 ' at w = 1/2 could be 
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found and because a Landau integrand singularity is a necessary condition^ 7 for an 
integral to be singular, it was conjectured that x^ 5 \ necessarily defined by analytic 
continuation from its series representation, would not be singular at w = l/sjjj. We 
have now, in section 4, shown that x^ 1S singular at w = 1/2 and so clearly have 
identified all the singularities of x^ by an analysis of the ODE it satisfies. The 
w = 1/2 singularity in x^ proves the existence of an associated Landau integrand 
singularity, but we have made no progress in identifying this Landau singularity. And 
so our goal of unifying the Landau integrand analysis with the ODE approach to the 
behaviour of the general x continues to elude us. 
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Appendix A. The exact arithmetic i^ 2 Ight ^ = L\ ■ Lu 

It was reported in [9] that = L\ ■ Lu with a solution S\2 that satisfies 

L\\ (S12) = P, L\(P) = where P is a polynomial. The order eleven linear 
differential operator Lu has a direct-sum decomposition 

L u = (Z 2 ■ Ni) ®V 2 ® (F 3 ■ F2 ■ LI) . (A.l) 

where Z2 is a second order operator also occurring in the factorization of the linear 
differential operator associated with x}- 3 ': V2 is a second order operator equivalent to 
the second order operator associated with x' 2 ' an d F3 and F2 are known in exact 
arithmetic [9]. The exact Lu was obtained in [9] from (|A.1[) and is giverjjj] in |10) . 
Our goal here is to determine the exact L^^ or equivalently the exact solution S12 
which gives the polynomial P = Lu(S 12) and then L\. 
One possible form for this solution is 

S12 = hx(w) ■ S 9 - 2w 3 - 6w 4 ~ 136 w 5 /3 ~ 3161 w 6 /24 

- 3073u> 7 /4 - 245147 w 8 / 12 - 334031 w 9 /30 

- 872442 u) 10 /35 - 19192553 u; n /140 + 0{w 13 ) (A.2) 

where 

S g = w 2 +4w 3 +30w 4 +120u; 5 + 690u; 6 + 2760w 7 + 14280 w 8 
+ 57120 w 9 + 279090 w 10 + 1116360 w 11 

+ 5261244 w 12 + 0(w 13 ) (A.3) 

§ Hwa and Teplitz (see page 6 in [7]) give 1141 and 1 151 as the original sources for the necessary 
requirement. 

I The abstract in [8] erroneously stated the absence of a singularity without caveat. In view of the 
findings of this paper the conclusion that \^ is singularity free at w 2 = 1/8 must now also be 
discounted. 

% Note that Lu in I10| is given in monic form. Here Lu is understood to be normalized to a sum 
of products of polynomials times powers of w ■ d/dw. 
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is annihilated by the order nine operator V 2 Ni (F 3 ■ P 2 • Lf). The rational 
coefficients in (|A.2[) are sufficiently small that the four modulo prime solutions of 
^ reported in [9] will yield S12 exactly to 0(w 19 ). This in turn means that 
P = Ln(5i2) modulo any prime is given exactly to 0(w 19 ). We now observe 
that the degree of P depends on the In representation and is 18 when either the 
minimum degree or minimum degree plus one representation of L\\ is used. This 
implies that, with this degree restriction on L\\ and without any further input from 
the equation Pis = Lu(Si 2 ) determines Pis modulo any prime exactly. The 
same equation can then be used in recursion mode to extend the known S12 modulo a 
prime series of 0(w 19 ) to any desired length. Clearly these steps can be repeated for 
as many primes as necessary to get the S12 coefficients in exact arithmetic by rational 
reconstruction. Then, once S12 has been extended to 0(u> 137 ), one can switch back to 
the minimum order representation of L u to get P136 = IiiO?^), followed by Li and 
the minimum order L^F^ in exact arithmetic. Note that we make no attempt to 
reconstruct Pis in exact arithmetic. This polynomial is of no particular use and since 
it is generated from a non-minimal order Lu we expect its coefficients to be equally 
outrageously large. 

The polynomial P136 is proportional to 

5544 w 5 + 197765 w 6 + 419883469 w 7 /70 - 24564638026 w & /35 + ••• 

-2 171 • 134196985 w 135 /363 - 2 173 • 80400525 w 136 /847. (A.4) 

and the operator L\ that annihilates it is obviously proportional to Pi36^ — d d™ • 
Straightforward operator multiplication L\ ■ Lu yields an 4 r 2 ight) and this can be 
normalized to monic or non-monic form as desired. Note that Piss/uP becomes 
the apparent singularity part of the head polynomial of L^ r 2 lsht ' . The corresponding 
apparent polynomial P a pp of Lu divides out of the product L\ ■ L\\ as a common 
factor; this serves as a useful check of our algebra. The minimum order L^^ we 
have obtained can be found in [TO] . 

The fact that the series Sg multiplying the logarithm in (|A.2[) is annihilated by 
an operator of order M = 9 < 11 suggests that it might be possible to "push" either 
V2 or Z2 to the left of L\. And indeed we found that a U\, equivalent to L±, occurs at 
the left of the L 9 = V 2 Ni © (P 3 • F 2 ■ If) that annihilates Sg in (|A~3)l . Next, our aim 
is to see whether U± can be pushed further to the right. This amounts to considering 
the various factorizations of Lg; for instance the factorization where the equivalent 
of V2 occurs at the leftmost position and there exists a direct sum of this with Ui. 
There is such a direct sum showing that an equivalent to U\ occurs at the left of 
Ni © (P3 • P2 ■ The process continues until it is no longer possible to find a direct 
sum of our equivalent to L\ and any one of the right factors. This final factorization 
has a solution 

S 5 = ln(w) ■ w 2 ■ l + V(l-4w)/(l+4w) /(l-4w) 2 

+64 w 5 /9 + 200 w 6 /9 + 712 w 7 /3 + 10708 w 8 /!5 
+253888 w 9 /45 + 1771852 w 10 /105 + 12245672 w n /l05 
+331735612 w 12 /945 + 0(> 13 ) (A.5) 

which can replace S12 from (|A.2[) as one of the 12 linearly independent solutions of 
L^ 2 lsht ^. The operator that annihilates S5 is 

U s = £1 • (JVi (F 2 ■ LD) (A.6) 
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where if we define U4 = Ni © (F 2 ■ L\) to be given in our standard non-monic form 
then U 4 (S 5 ) oc P 19 w 2 /(l - Aw) with 



P19 = 9 + 36 w + 18 w 2 - 2064 w 3 + 4581 w 4 + 59584 w 5 - 143476 - 898464 w 7 
+124724 w 8 + 813120 w 9 + 9220240 w w + 55704896 w 11 + 65556224 w 12 
-253883392 w 13 - 406194176 w 14 + 1318182912 w 15 + 2053013504 w 16 
+368443392 w 17 - 454033408 w ls - 272629760 w 19 . (A.7) 

The operator that annihilates Kt(Ss 

u 



>5j IS 
d 

dw 



(2~Aw)P 19 + w(l-A 



dP 19 
dw 



(A.8) 



and again it is straightforward to combine this L\ with the other exact arithmetic 
operators in (|A.6[) to get U5. 

An incidental remark on the solution (|A.5[) is that this is one of the few 



examples of a (partial) analytic solution of L 



(right) 



Other simple solutions are 



w 2 /(l - Aw) 2 , w 2 /(l - 4w)/y/l - I6w 2 and w 2 (1 + 4 w) 2 F X (3/2, 3/2; 1; 16 w 2 ), 
w 2 (1 +4w) 2 Fi(3/2,3/2;3;l- 16w 2 ) associated with Lf, Ni and V 2 , respectively. 
A more complicated case is that of Z 2 for which we refea the reader to [9] . 

In conclusion, every solution of L^" sht ^ can now be obtained as a linear 
combination of solutions of operators of order M < 6. 



Appendix B. a * the ferromagnetic point 

To find $( 5 ) in the vicinity of the point y = where 2y = 1 —Aw requires that we first 
find a solution to the inhomogeneous equation ([2]) near y = 0. To this end we develop 
the elliptic integrals in Q as series in y after which finding a particular integral by 
series recursion is straightforward. What remains is then the matching of the <fr( 5 ) 
series about w = to this particular integral plus the 24 homogeneous solution series 
Si satisfying £24(1%) = 0. This too is a straightforward exercise although clearly 
one must pay attention to potential errors arising both from numerical round-off and 
series length truncation. A useful remark in this regard is that it is very advantageous 
to work with <I>( 5 ' series in s rather than w — s/2/(l + s 2 ). This arises because of the 
non-linear relationship 1 — Aw w (1 — s) 2 /2 at the ferromagnetic critical point. 

To understand the effect of the choice of series variable we note that in general 
the error resulting from truncation of a series in x at terms when evaluated at 
x m scales roughly as (x m /x s ) N where x s is the distance to the nearest singularity 
froirfj] x = 0. As an illustration, from $( 5 ) series of N = N w terms and a 
matching point at w m = (1 — 2y. m )/4, the series truncation error is w (Aw m ) N 
= (1 — 2y m ) N in exp(—2N w y m ) for small y m . Now suppose the solution series 
Si have been evaluated to N = N y terms. The solutions with the smallest radius 
of convergenc^J, y s = 1/2 — 2w s m 0.183, determine the truncation error which is 

1 Note a misprint in Appendix C in [9], in eqn (C7). The right-hand side should in fact be replaced 
by 8 2 (1 - 2w) 2 (l + 2w) 2 /(l - w) 2 /(l - Aw)' 2 . 

t We assume for purposes of our qualitative discussion here that all functions are approximately of 
unit magnitude. 

\ This w s is the same as that in the discussion leading to and following ||9I10( . 
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~ {y m /0-^83) N ~ exp(— N y ln(0.183/y m )). The optimal choice of matching point is 
where the errors are roughly equal, i.e. 

2N w -y m = N y • ln(0.183/y m ). (B.l) 

With N w = 8000 and N y = 800, condition^] ([BT]l yields y m = 0.0577 with a 
resulting truncation error w 10~ 400 . This estimate of 400 digit accuracy is overly 
optimistic since the <£>( 5 ) matching requires the evaluation of 23 derivatives which 
invalidates our assumption of functions all of unit magnitude. 

Consider now what happens if the matching is done using $^ 5 ^ as a series in 
s of N = N s = N w terms. The truncation error is ~ {s m ) N ~ (1 — 2 y/y m ) N 
~ exp(— 2 N w y/ym) for small y rn . The match point choice that replaces (|B. 1[) is 

2N W ■ ^ = N y ■ ln(0.183/2/ m ), (B.2) 

and with the same series lengths as above, y m = 0.01535 with resulting truncation 
error w 10~ 860 . Again this must be overly optimistic but it does illustrate the 
dramatic improvement achieved with no extra computational cost. 

Our results for the behaviour of $' 5 - ) at the ferromagnetic point are as follows: 
We use throughout the definition 2y = 1 — 4w as above. Then, to 0(y 10 ), the 
5-particle contribution is 

120tt 4 -$( 5 )= (x (1) - 60x (3) + 120x (5) W 4 = 
ln 4 (y/4) • (-5/128 - 25y/512 - 75y 2 /512 + 6455?/ 3 /8192 

+42305y 4 /16384 + 49935y 5 /8192 + 444955y 6 /32768 

+64409075 2/ 7 /2097152 + 306977235 y 8 /4194304 

+3131527805y 9 / 1 6777216) 
+ ln 3 (y/4) • (-235/192 - 1043y/512 - 559 y 2 /256 

+1007843 y 3 /1032192 + 265627475 y 4 /22708224 

+21838013489 y 5 /590413824 + 57270165499 ?/ 6 /590413824 

+315303684751873 y 7 / 12 84740481024 

+31001275613111851 ?/ 8 /48820138278912 

+336259319399213305 y 9 / 1 95280553115648) 
+ ln 2 (y/4) • (-1225/192 + 306283 y/15360 - 859553 y 2 / 215040 

-20921323001 y 3 /433520640 - 10704883015027 ?/ 4 /104911994880 

-12546069917407919 y 5 /70920508538880 

- 1021232328273251 y 6 /3546025426944 

-1219521777522960411013 y 7 / 2 623491451870248960 

-1428705274437356046530557 y 8 /1894160828250319749120 

-3094013953756589876579173 y 9 /2525547771000426332160) 

+ ln(y/4) • • y k j + ]T C (k) ■ y k (B.3) 

\fc=0 / k=-l 

where the constant arrays C\ and Co are only known as floating point values. 
Truncated values are 

Ci(0) = -24.57942277608500580980766691884235675672213367715281 • ■ ■ 

1 We have reached such values for exact rational coefficient series without too much computing effort. 
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c. 


'1 1 




1UU.00ZZ0UZU1J.Z1 you i Oc/OllOO ( O 1 OUOyU4:Z1Z f IIIUOIO ± 1 uuzuo 


C\ 
L'l 


+ J 




1 £901 999^00^00987037^/1981 37031 Q1 7^0/170(^^^700/190/1/1 1 
loo.uzy izzzouyooyuzo / yo / u^zoio < uoiyi / oy^ / yuoo / yy^zy^^ii • • • 


Ci 


'3) 




96.268269775450223433437918344019644004451927310385706 • • • 


Ci 


;4) 




-93.06087224694711987327189827852646701359704351289178 • • 


Ci 


:s) 




-520.0604636206711421998793353361815107994358750212855 • • 


Ci 


;o) 




-1469.482403857530493377284131808565295376128516174659 • • 


Ci 


:t) 




-3700.403260075047090620399746433870069155277603230413 ■ • 


Ci 


18) 




-9303.076094072409117026994688305643439126188417678639 • • 


Ci 


'9) 




-24403.45115956242143888235859126942238116466146718368 • • 


^0 




— 93 1 fi4^fi1 9033fifi71 91 1 7448^.48000^08800004398948^1 0fi70fi01 








— 117 ^^1 74nfi93n09343n801 0^781 40^,81 4973001 R3833071 fi91 . . 


<~-0 


:i) 




9^R 41 1 ^1 1 40Q40fi939^7Q981 OH31 49Q34fi1 1 01 8878000^93401 8 . . . 
zou.iiioii^yy^yozozo/yzoiuuoi^zyo^uiiyioo/ouyuouzo'iyio 


<~-0 


'2) 




3^0 01 ^.8^.00^7001 01 ^9091 08301 1 fifil Q800RQ1 81 90fi3077fi1 1 000 . . . 

oou.yioooyuo /yuiuiozyziyooyiiuuiyouuoyioizyuoy ( iui±yyy • • • 




[3) 




903 A 780301 K7K81 8fi30Rfi9393fi3f;0fifi8073Q1 ^.1 4 34 08^01 4 8^7/1 1 
Zyo. 4 I oyoy ID 1 DolouoUDDZoZoDoDUuDoU (oy 1014041100111400(41 • • • 


Co 


;4) 




109.38934363323783704769948377240709436366681466990651 ■ • • 


Co 


:s) 




-373.9288055773641613398692905165600188023113404992740 • • 


Co 


;o) 




-1561.022604360024268401667883538021798990658115951832 ■ ■ 


Co 


:t) 




-4574.248865947314517242004568453306010341709549402520 • • 


Co 


;») 




-12595.46119191674598711644766973884064574672223482538 • • 


Co 


'9) 




-35224.55961370779008532569450932078801945477728287548 • • 



(B.4) 

while values to about 800 digits can be found in [ID]. The 3-particle function, to 
0(y m ), is 

-6^ 2 -$( 3 > = -6x (3) )-tt 2 = 

ln 2 (y/4) ■ (3/32 - 3y/128 + 9y 2 /128 + 87y 3 /2048 

-39y 4 /4096 - 423y 5 /4096 - 555//2048 - 301647 y 7 /524288 

-1185363 y 8 /!048576- 8996817 //4194304) 
+ ln(y/4) • (23/32 + 227y/1280 - 2047y 2 /4480 - 88949 y 3 /122880 

-20562503 y 4 /!8923520 - 810591833 y 5 /492011520 

-316471567 y 6 / 12 3002880 - 4445362809179 y 7 /1070617067520 

-3085016829083 y 8 /447070863360 

-273049228448281 y 9 /23247684894720) 
+ (41/96 - 21169^/38400 - 225583 y 2 /3763200 + 2621231 y 3 / 154828800 

+15533081173 y 4 /262279987200 

+151705656477979 y 5 /1595711442124800 

+30692556260057 y 6 /265951907020800 

+1747927197871890533 y 7 /19676185889026867200 

-41073304786882831381 y 8 /655671055932802990080 

-776816275820131600824097 y 9 /!534270270882758996787200) 
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+ (9V (1) (!/3)/8 - 3/2 - 3tt 2 /4) • (l/y-3/2 - 5y/24 - 3y 2 /8 
-1801 y 3 /3456- 1649 y 4 /2304 - 187999 2/ 5 /!86624 
-45617 y 6 /3H04 - 17592665?/ 7 /7962624 

-164030851 y 8 /47775744 - 9425604977y 9 /m9926784) (B.5) 

where ij)^ l '{x) = ifi'(x) is the polygamma function (ip(x) — T'(x)/r(x)). The 
coefficient of 1 jy in (|B.5[) was determined by Tracy [16] while the higher order terms 
are recursively generated from the ODE satisfied [T3] by $( 3 ). Finally, 

= x (1) = V(%) - 1/2. (B.6) 

When these results are combined and expressed in terms of r = (1/s — s)/2, i.e. 
2 y = t 2 /(1 + t 2 + \/T+ t 2 ), one confirms the low order result for x^ m equation 
(48) of [8]. We also confirm the 500 digit expression for the x^ amplitude D§ given 
by Bailey et al [TT]. The explicit connection is 

D 5 = 52tt 4 /5 - 12tt 2 i/' (1) ( 1 /3) + 16tt 2 + 16C (-1)/15 (B.7) 

with C (-l) given in (|B.4|) . 

The 5-particle contribution (|B~3l) to 0(y w ) can be extended to arbitrary order by 
a purely local analysis of the ODE @. This makes (|B.3|) particularly useful for those 
problems in which one wants to analytically continue <i>( 5 ) beyond the ferromagnetic 
singularity. Of the constants in (|B.4|) . only the fifteen Ci(0..5) and Cq(— 1..7) depend 
directly on solution matching. All higher order coefficients are fixed by recursion 
relations with rational coefficients. The rationals are quite large, even in the simplest 
example which reads 

28235138014521201151215868809287621440- C x (6) = - — 

-37737948182868464640476571846064977272634973/1168128 • d(0) 
-2562052253560544220030628882123379039659169/86528 • d(l) 
+8187149177658628591337051452109148482747/312 • Ci(2) 
+158617170239026490613752192866932708861763/1872- Ci (3) 
-140695455847623820725386797097484481682383/702 • C X {A) 
+ 133960461264493957765301921009624203227 • C x (5), (B.8) 

with 

m = 15256931516199856571741494453893800774223260951365779571949, 
di = 31365304106405068800. 

Floating point recursion in y is unstable to the extent that relative errors increase 
by about factor 40 per order in the y series. If one wants, say, 1500 terms, one must 
start with 2400 digits more than the desired final accuracy. One must also use directly 
as starting values only the coefficients Ci(0..5) and Co(— 1..7) from (|B.4|) . All other 
coefficients are to be obtained to the necessary higher accuracy directly from the ODE 
^ locally around y = 0, now treating the starting values Ci(0..5) and Co(— 1..7) as 
if exact. 

The problem of numerical instability in series generation around rational points 
can be avoided by generating all series in exact arithmetic. In that case the series 
coefficients are in general rational and our observation, admittedly rather limited, 
is that the growth of the denominators is at most exponential. An example of 
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Figure Bl. The number of digits D in the denominators of the rational 
coefficients cjv in series y p In 9 (y / 4) ^2c^y N that arise in ODE solutions around 
2y = 1 - Aw = 0. 



denominator growth is shown in Fig. IB1I for all series in y that we have generated 
for doing the ferromagnetic point series matching described above. A growth rate 
faster than exponential would make exact arithmetic series generation impractical 
so it is of some importance to know whether the observed exponential growth is a 
general feature. It is reminiscent of G-series |17) but we are unaware of any theorem 
that guarantees this behaviour for expansions of integrals with algebraic integrands 
about singular points. 

Appendix C. Right division operators for irrational singular points 

We wish to determine a right division of a minimum order Ln(w) from solutions 
identified with the non-rational zeros of a factor h(w) of the head polynomial of 
Ljy(w). Examples of h{w) treated in the text are 1 + 3w +4w 2 and the case 3 
and 4 polynomials from Appendix C in \8\ . If such a polynomial has an "accidental" 
modulo prime p factorization h(w) = (w — w p ) h p {w), we can still generate series 
solutions about x = where x — w — w p . Furthermore, suppose the solution S^cc) 
has been tested and is found to be annihilated by a minimum order M < N operator. 
If the code used for the testing is that described in section 3 of [5] the result will be 
of the form 



where the head polynomial fo{x) is normalized such that f (x — 0) = 1. Note that 
Wp is not the modulo p representation of a root of h(w) and both it and x have no 
meaning independent of the specific modulo p factorization h(w) = (w — w p ) h p (w). 




(C.l) 
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Consequently there does not exist any exact ll Lj^{xf to be reconstructed from 
multiple L M (x) generated using different primes. Replacing x by w — w p in (|C 1[) 
generates an Lm{w) but only when this is correctly normalized to remove all 
dependence on the "accidental" w p can it be used for a reconstruction of an exact 
Lm(w). The steps for doing this are as follows. 

To begin, note that (|C.1|) can be expressed as a sum of polynomials times powers of 
d/dx. The head polynomial in this representation is clearly x fo(x). The modulo p 
factorization of this will contain as factors x M [h p (x + w p )] K , typically with K < M. 
These factors can be written as x M ~ K [xh p (x + w p )] K — x M ~ K [h(x + w p )] K = 
(w ~ w p ) M ~ K {h(w)] K . Dependence on the "accidental" w p appears only in the factor 
x m-k an( j can be divided out. So too can the w = normalization constant 
C = x K fo(x) evaluated at x = —w p . The result is the right division (modulo p) 
operator which can be cast into the same form as (IC.1|) . namely 

Lm{x) 



L m {w) = , f '_y modulo p 



d 



= £^)>-^J (C2) 

71 = 

where each g n (w) is a polynomial and, by construction, go(w = 0) = 1. The exact 
Lm(w) can now be reconstructed from different prime based versions of (|C.2[) . 
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